Ultimate Fate of Constrained Voters 
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We determine the ultimate fate of individual opinions in a socially-interacting population of 
leftists, centrists, and rightists. In an elemental interaction between agents, a centrist and a leftist 
can become both centrists or both become leftists with equal rates (and similarly for a centrist and a 
rightist). However leftists and rightists do not interact. This interaction step between pairs of agents 
is applied repeatedly until the system can no longer evolve. In the mean-field limit, we determine 
the exact probability that the system reaches consensus (either leftist, rightist, or centrist) or a 
frozen mixture of leftists and rightists as a function of the initial composition of the population. We 
also determine the mean time until the final state is reached. Some implications of our results for 
the ultimate fate in a limit of the Axehod model are discussed. 
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I. INTRODUCTION 

A basic issue in social dynamics is to understand how 
opinion diversity arises when interactions between indi- 
viduals are primarily "ferromagnetic" in character. Many 
kinetic spin models have been proposed to address this 
general question 0, IH, 0, || ■ An important example of 
this genre is the appealingly simple Axehod model |(| Q , 
which accounts for the formation of distinct cultural do- 
mains within a population. In the Axehod model, each 
individual is endowed with a set of features (such as polit- 
ical leaning, music preference, choice of newspaper, etc.), 
with a fixed number of choices for each feature. Evolu- 
tion occurs by the following voter-model-like update step 
[j| . A random individual is picked and this person selects 
an interaction partner (anybody in the mean- field limit, 
and a nearest-neighbor for finite spatial dimension). For 
this pair of agents a feature is randomly selected. If these 
agents have the same state for this feature, then another 
feature is picked and the initial person adopts the state 
of this new feature from the interaction partner. This 
dynamics mimics the feature that individuals who share 
similar sentiments on lifestyle issues can have a meaning- 
ful interaction in which one will convince the other of a 
preference on an issue where disagreement exists. 

Depending on the number of traits and the number of 
states per trait, a population may evolve to global con- 
sensus or it may break up into distinct cultural domains, 
in which individuals in different domains do not have any 
common traits 0, Q • This diversity is perhaps the most 
striking feature of the Axelrod model. An even simpler 
example with a relat ed p henomenology is the bounded 
compromise model yj, [I(J, . Here, each individual pos- 
sesses a single real-valued opinion that evolves by com- 
promise. In an update step, two interacting individuals 
average their opinions if their opinion difference is within 



a pre-set threshold. However, if this difference is greater 
than the threshold, there is no interaction. These steps 
are repeated until the system reaches a final state. For 
a sufficiently large threshold the final state is consensus, 
while for a smaller threshold the system breaks up into 
distinct opinion clusters; these are analogous to the cul- 
tural domains of the Axelrod model. 

In spite of the simplicity of these models, most of our 
understanding stems from simulation results. As a first 
step toward analytic insight, a discrete three-state ver- 
sion of the bounded compromise model was recently in- 
troduced This is perhaps the simplest opinion dy- 
namics model that includes the competing features of 
consensus and incompatibility. This model was found to 
exhibit a variety of anomalous features in low dimensions, 
including slow non-universal kinetics and power-law spa- 
tial organization of single-opinion domains. In this work, 
we focus on the ultimate fate of this system in the mean- 
field limit. We determine the exact probability that the 
final state of the system is either consensus or a frozen 
mixture of leftists and rightists as a function of the ini- 
tial composition of the system. We also compute the time 
required for the system to reach its ultimate state. 

In the next section, we define our model, outline its 
basic properties, and determine the ultimate fate of the 
system in terms of an equivalent first-passage process. 
The solution to this problem is used to obtain the prob- 
abilities of reaching either a frozen final state (Sec. Ill) 
or consensus (Sec. IV). In Sec. V, we compute the mean 
time until the final state is reached. In the discussion 
section, we generalize to non-symmetric interactions and 
also show how our results can be adapted to determine 
the ultimate fate of a version of the Axelrod model. Cal- 
culational details are presented in the appendices. 



II. THE MODEL 
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We consider a population of N individuals, or agents, 
that are located at the nodes of a graph. Each agent 
can be in one of three opinion states: leftist, centrist, or 
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rightist. As shown in Fig.Q] we represent these states as 
— , 0, and +, respectively. In a single microscopic event an 
agent is selected at random. We consider the mean-field 
limit in which the neighbor of an agent can be anyone else 
in the system. If the two agents have the same opinion, 
nothing happens. If one is a centrist and the other is 
an extremist (or vice versa), the initial agent adopts the 
opinion of its neighbor — that is, each individual can be 
viewed as having zero self-confidence and merely adopts 
the state of a compatible neighbor; this kinetic step is 
the same as in the classical voter model ||. However, 
if the two agents are extremists of opposite persuasions, 
they are incompatible and do not influence each other. 
As a result of this incompatibility constraint, the final 
state of the system can be either consensus of any of the 
three species, or a frozen mixture of leftist and rightist 
extremists, with no centrists. 
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FIG. 1: Update events for different pair states. 

In a mean-field system with A_ leftists , A + right- 
ists, and A centrists (with A+ + A_ + A = A), the 
probability of selecting a pair ij = +,—,0) equals 
2NiNj/[N(N - I)], where N{ is the number of agents 
of type i. Then the elemental update steps and their 
respective probabilities are: 



(A_,A )^ (A_±I,A Tl) prob. p x 
(N+,N ) -> (A + ±l,A oT l) prob. Py 



A_A 
A(A- 1) 

A+A 
N(N - I) : 



(1) 



while the probability for no change is 1 — 2(p x + p y ). 

Eventually the system reaches a static final state that 
is either consensus of one of the three species or a frozen 
mixture of leftists and rightists. Monte Carlo simulations 
showed that the nature of the final state has a non-trivial 
dependence on the initial densities of the species [l^ . 
We now analytically determine this final state probability 
as a function of the initial population composition by 
solving an equivalent first-passage problem. As shown 
in Fig. |2 the state of the system corresponds to a point 
in the space of densities x = A_/A, y = A + /A, and 
z = Aq/A. Since x + y + z = 1, this constraint restricts 
the system to the triangle ABD shown in the figure. 

When two agents interact, the state of the system may 
change and we can view this change as a step of a corre- 
sponding random walk on the triangle ABD, with single 
step hopping probabilities given in Eq. When the 
walk reaches one of the fixed points A, B, or D (con- 




FIG. 2: The density triangle x + y + z = 1. The heavy dots 
denote consensus states and the heavy line denotes the frozen 
final states where no centrists remain. Typical random walk 
trajectories are shown, along with their projections (dashed) 
on the z = plane. When a random walk reaches the lines 
x = or y = 0, the random walk subsequently must remain 
on this line until consensus is reached. The corresponding 
hopping probabilities on the z = plane are also indicated. 



sensus), or any point on the fixed line AB (frozen ex- 
tremist mixture), the system stops evolving. This set 
represents absorbing boundaries for the effective random 
walk. Thus to find the probability of reaching a given 
final state, we compute the first-passage probability for 
the random walk to hit a given absorbing boundary. 

We may simplify the problem by using the fact that 
only two of the densities (x,y,z) are independent. We 
thus choose x and y as the independent variables by 
projecting the effective random walk trajectory onto the 
z = plane (Fig. According to Eq. this two- 
dimensional random walk jumps to its nearest neighbors 
in the x and y directions with respective probabilities 
p x = xz/(l — d) ~ xz and p y — yz/(l — 6) w yz, and 
stays in the same site with probability 1 — 2(p x + p y ). 
Here 8 — 1/A, and we consider the limit A — > oo. 



III. FROZEN FINAL STATE 

We now determine the probability {x,y) for the 

system to reach a frozen state when the initial densities 
are p-(i = 0) = x and p + (t = 0) = y. This coincides 
with the first-passage probability for the equivalent ran- 
dom walk to hit the line x + y = 1 when it starts at a 
general point (x, y) in the interior of the triangle. There 
are also absorbing boundaries on the sides x = and 
y = 0, where the probability of reaching the frozen state 

is zero, and P_| (x, y) = 1 on the line segment x + y = 1. 

For notational simplicity, we write F(x,y) = (x,y). 

In the equivalent random walk process, this first-passage 
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probability obeys the recursion 01 

F(x,y) = Px [F{x-S,y) + F{x + S,y)} 
+ p y {F(x,y-S) + F(x,y + S)} 
+ [l-2( Px +p y )}F{x,y). 



(2) 



That is, the first-passage probability F{x, y) equals the 
probability of taking a step in some direction (the factor 
Pi) times the first-passage probability from this target 
site to the boundary. This product is then summed over 
all possible target sites after one step of the walk. 

In the continuum limit (N — > oo) we expand this equa- 
tion to second order in 5 and obtain 



(3) 



x d 2 F(x,y) | d 2 F(x,y) _ Q 
dx 2 dy 2 

supplemented with the boundary conditions 

F(x,0) = 
F(0,y) =0 
F(x,l-x) = 1. 

The solution to Eq. is not straightforward because 
of the mixed boundary conditions. Instead of attacking 
the problem directly, we transform to the coordinates 
u = i/x, v = ./y to map the triangle to a quarter-circle 
of unit radius and then apply separation of variables in 
this geometry. 

The transformed differential equation is 



d 2 F(u,v) d 2 F{u,v) 



du 2 



dv 2 

ldF(u,v) ldF(u,v) 
u 



du 



dv 



0. (4) 



Because of the circular symmetry of the boundary con- 
ditions in this reference frame, it is now convenient to 
use the polar coordinates (p, 9) , where u = p cos 9 and 
v = psind. This transforms Eq. Q to 



d 2 F{p,( 
dp 2 



1 dF(p, i 



P dp 



+ — j (tan6> — cot ( 



I d 2 F{p,l 

ae 2 

dF{p,i 



09 



(5) 



We seek a product solution F(p, 9) = R(p)Q(9). Substi- 
tuting this into Eq. © leads to the separated equations 



d 2 R 
dp 2 

d 2 S 



1 dR m 2 

: 7rR = 



p dp 



(6) 

+ (tan6»-cot6»)^-+m 2 e = (7) 



dQ 



d9 



where m is the separation constant. 

Eq. 10 is equidimensional and thus has the general 
power-law form 



R(p) = A +P 1+VI +^ + A-p 1 -^ 1 ^ 



(8) 



where A± are constants. To solve Eq. Q, we proceed by 
eliminating the first derivative to give a Schrodinger-like 
equation. Thus we define <d(9) = f{9) u(9) and find the 
function f(9) that eliminates this first derivative term. 
Substituting 6 = / u in Eq. @ gives 

fu" + [2/' + (tan 9 - cot 9) /] v! 

+ [/" + (tan 9 - cot 9) f + m 2 f] u = 0, (9) 

and the coefficient of u' is zero if / satisfies 

2/' + (tan0-cot0)/ = O, 

whose solution is 



f(9) cx \/2sin0cos0 cx Vsin 29. 



(10) 



Substituting this expression for / in Eq. @, we obtain 



cPu 
dtP 



(1 + m 2 )-- 



sin 



0. (11) 



Details of the solution to this Schrodinger equation are 
given in Appendix A. The final result is 

F(p,9) = R{p)f(9)u{9) 

= E ^-^|p 2(n+1) S in20^(cos20), 

(12) 

where is the associated Legendre function. We trans- 
form back to the original xy coordinates through 



p cos 9 = 



sin 29 



x 
xy 



p sin 9 = 

x-y 



cos26> 



y 



x ■ 



y 



Finally, identifying F with P^ , the solution in the orig- 
inal Cartesian coordinates is 



x-y 
x + y 



(13) 

This gives the probability of reaching a final frozen 
state when the population has initial densities x and y of 
leftists and rightists, respectively. Notice that this prob- 
ability is symmetric, P^ (x,y) — P + _(y,x), because P^ 

is an even function for n odd and an odd function for n 
even. This reflects the obvious physical symmetry that 

the probability P^ is invariant under the interchange of 

leftists and rightists. 

For equal initial densities of leftists and rightists, 
namely x = y = (1 — z)/2, Eq. i|13|) becomes a function 
of the initial density of centrists z only 



n odd 



i^±n ( i-^+ipi ( o). 

n(n + 1) 
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As shown in Appendix B, this can be simplified to the 
closed-form expression 



(14) 



with < z < 1. This solution is shown in Fig. [21 along 
with Monte-Carlo simulation results. 
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FIG. 3: Probability of reaching the frozen final state, P+-, 
and extremist consensus, P+ , as a function of the initial den- 
sity of centrists z for the case of equal initial densities of 
leftists and rightists, y/x — 1 (top) and for the case y/x = 9 
(bottom). 



IV. CONSENSUS FINAL STATES 

In addition to reaching a frozen final state, the system 
can also reach global consensus — either leftist, right- 
ist, or centrist. Again, the probabilities for these latter 
three events depends on the initial composition of the 
system. We thus define P-(x,y), P + (x,y), and Po(x,y), 




FIG. 4: Probability distribution that an initially symmetric 
system (x = y) has magnetization m in the frozen final state, 
for the cases z = 0.04, 0.1, and 0.2 (progressively broadening 
curves) . 



as the respective probabilities to reach a leftist, rightist, 
or centrist consensus as a function of the initial densities. 

The probability for centrist consensus can be obtained 
by elementary means. We merely map the 3-state con- 
strained system onto the 2-state voter model by consider- 
ing the leftist and rightist opinions as comprising a single 
extremist state, while the centrist opinion state main- 
tains its identity. The dynamics of this system is exactly 
that of the classical 2-state voter model. Since the over- 
all magnetization of the voter model is conserved Q , the 
probability that the system reaches centrist consensus, 
for a given x and y is simply 



Po(x,y) = z = 1 - (x + y). 



(15) 



To find the probability for leftist and rightist consen- 
sus, we use normalization of the final state probability 

P+(x,y) + P-(x,y) + P+-(x,y) + P (x,y) = 1 (16) 

and conservation of the global magnetization to write 



y-x = P + (x,y) - P-(x,y) 



(1 - 2x')F(x'\x,y)dx'. 



(17) 



Here F(x'\x, y) is the probability of ending in a specific 
frozen state with density x' of — spins and 1 — x' of + 
spins as a function of the initial densities x and y. This 
function therefore satisfies the normalization condition 

P_i (x,y) — J Q F(x'\x, y)dx' . Additionally the integral 

in Eq. (|17l) is the final magnetization in the frozen state. 
From the exact solution for F(x'\x, y) given in Eq. I|A6|I . 
we thereby obtain as a byproduct the magnetization dis- 
tribution in the final frozen state (Fig.^J. As expected, 
for a small initial density of centrists z, there is little 
evolution before the final state is reached and the mag- 
netization distribution is narrow. 
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Using Eq. (|15fl and the normalization condition for 
F(x'\x,y), we recast Eqs. ljTtj|) and lfT7|l as 

P+(x,y) + P^{x,y) + P + _(x,y) = x + y, (18) 

and 

P+(x,y) - P-(x,y)+P+-(x,y) 

-2 x'F(x'\x,y)dx' = y-x. (19) 
Jo 

Subtracting these equations, we obtain 

P-(x,y) =x- [ xF(x'\x,y)dx' (20) 
Jo 

Now the first-passage probability F(x'\x, y) to a frozen 
state with a specified density x' of — spins obeys the 

same differential equation as P_| [Eq. (J3J], but with the 

boundary conditions 

F(x'\x,0) = 
F(x'\0,y) = 
F(x'\x,l-x) = S(x-x') (21) 

The last condition states that the first-passage probabil- 
ity to the point (x' ,1 — x') on the boundary x + y = 1 is 
zero unless the effective random walk starts at (x' , 1—x'). 
The first-passage probability F(a/|:r, y) has the general 
form given in Eq. (|A4|) . but with the coefficients c„ now 
determined from the boundary conditions in Eq. (|2Tjl . 

From the solution for F(x'\x, y) given in Appendix A, 
the probability of — consensus as a function of the initial 
densities x and y is 

P {x,y)=x-jr ^±±> (a + y) n P l n ( X —A 
n(n + 1) \x 4- y ) 

71=1 V 

(22) 

The probability of + consensus can be obtained by using 
the fact that P+(x, y) = P-(y, x). 

For the special case ofx = y = (l — z)/2, the proba- 
bilities P+{z) = P-(z) can be obtained either by setting 
x = y in Eq. 12 I'D and summing the series, or, more sim- 
ply, by using Eq. (|14fl for P^ (z) and the probability 

conservation equation 2P + (z) + P^ (z) + z — 1. By ei- 
ther approach, we find the following closed expression for 
the probability of extremist consensus as a function of z 

P+(z) = P.(z) = l( l -( l ~ ^L-z\. (23) 

This result is also shown in Fig. [21 

From our results for Pj, (x,y), P + (x,y), P- (x,y), and 

Po(x, y), we may define a "phase diagram" shown in fig- 
ure El Each region in the figure corresponds to the por- 
tion of the composition triangle where the probability of 
ultimately ending up in the labeled state is greater than 
all other first-passage probabilities. The least likely out- 
come is the achievement of extremist consensus while the 
most likely result is to get stuck in the frozen mixed state. 



D 




FIG. 5: Exact phase diagram in the composition triangle 
ABD of Fig. |21 The triple points, where the final state prob- 
abilities of three of the four phases are equal, are located at 
(x,y) « (0.0507,0.6185) and (x,y) w (0.6185,0.0507). 



V. MEAN EXIT TIME 

In addition to the probability of reaching a particular 
final state, we also study the mean time until the final 
state is reached as a function of the initial composition of 
the system. The simplest such quantity is the uncondi- 
tional mean time t{x 1 y) to reach any of the four possible 
final states — extremist consensus (+ or — ), centrist con- 
sensus, and mixed frozen state, as a function of the initial 
densities x and y. This first-passage time can again be 
obtained trivially by considering an equivalent two-state 
system where we lump leftists and rightists into a sin- 
gle extremist state. The system stops evolving when the 
density of centrists reaches z = or z = 1, with the for- 
mer corresponding either to extremist consensus or to a 
frozen mixed state. 

To find this first-passage time to reach the final state, 
note that in a single event the effective one-dimensional 
random walk that corresponds to the state of the system 
can either jump to one of its two nearest neighbors with 
probability p z or stay at the same site with probability 
1 - 2p z , where p z = p x + p y = z(x + y) = z(l - z). The 
time interval for each event is dt = 1/N, corresponding 
to each person being selected once on average every N 
update steps. Then the mean time to reach the final 
states z = 0orz = lasa function of the initial density 
z obeys the recursion 

t(z) = p z [t(z + 6) + dt] + p z [t(z - 6) + dt] (24) 
+ {1 - 2p z )[t(z) + dt]. (25) 

This formula has a similar form and a similar explanation 
as the equation for the first-passage probability [Eq. (0)]. 
Starting from z, the mean time to reach the final state 
equals the probability of taking a single step (the factors 
p z and 1 — 2p z ) multiplied by the time needed to reach the 
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boundaries via this intermediate site. This path-specific 
time is just the mean first-passage time from the inter- 
mediate site plus the time dt for the initial step. 
In the large- N limit this recursion reduces to 



dH(z) 
dz 2 



1 



(26) 



where the diffusion coefficient is D = p z S 2 /dt. Since dt = 
S = 1/N, we have D = z(l — z)/N. Eq. is subject to 
the boundary conditions i(0) = t(l) — 0, corresponding 
to immediate absorption if the random walk starts at the 
boundary. The solution to this equation is 



t( z ) =-N[z Inz + (1 - z) ln(l - z)] 



(27) 



This is simply the mean consensus time of the 2-state 
voter model in the mean-field limit, in which the initial 
density of the two species are z and 1 — z. 
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FIG. 6: Unconditional and conditional exit times as a function 
of 2 for the cases y/x = 1 (top) and y/x = 2 (bottom). The 
curves for t and to are based on exact formulae, while the 
other are based on simulations. The latter results are based 
on 10 realizations for each initial state. 

It is more interesting to consider the conditional first- 
passage time to a specific final state as a function of the 
initial condition. For example, consider the conditional 



time to(z), the mean time to reach centrist consensus 
when the initial centrist density is z. This is the mean 
time for the equivalent random walk to hit the point 
z = 1, without hitting z = 0. In the large N limit, 
this conditional exit time obeys the differential equation 
(see Sec. 1.6 in 13]) 



D- 



[P (z) t Q {z)] = -P (z). 



(28) 



Integrating Eq. (|28|) . using Pq(z) = z, and the boundary 
conditions Pq(z) to(z) = at z = and at z = 1, we 
obtain Pq(z) t (z) — —N (1 — z) ln(l — z), from which 
the conditional time to reach centrist consensus is 



1 - z 

t Q (z) = -N ln(l - z). 



(29) 



Similarly, the mean time to reach all other final states of 
the system is 



Mother (z) = -N 



1 - z 



In z. 



(30) 



In principle, we can find the conditional times to reach 
the extremist consensus states and the frozen final state 
as a function of the initial condition. This involves solv- 
ing the two-dimensional analogue of Eq. (|28|l in the den- 
sity triangle, subject to the appropriate boundary con- 
ditions. Because of the tedious nature of the calcula- 
tion, we have instead resorted to numerical simulations 
to compute these conditional exit times. These results 
are shown in Fig. El While all exit times are of the order 
of N, it is worth noting that f_|_ and i_ are typically the 
longest times. This stems from the fact that reaching 
extremist consensus is a two-stage process. First, the ex- 
tremists of the opposite persuasion must be eliminated, 
and then there is a subsequent first-passage process in 
which the centrists are also eliminated. 



VI. DISCUSSION 

We determined basic properties of the final state in a 
simple opinion dynamics model that consists of a pop- 
ulation of leftists, rightists, and centrists. Centrists in- 
teract freely with extremists (either leftist or rightist), 
so that one agent adopts the opinion of its interaction 
partner, while extremists of opposite persuasions do not 
interact. These competing tendencies lead either to ulti- 
mate consensus or to a frozen final state that consists of 
non-interacting leftists and rightists. While this model is 
clearly an oversimplification of opinion evolution in a real 
population, it provides a minimalist description for how 
consensus or distinct cultural domains can be achieved. 

Our calculations are based on the mean-field limit in 
which each individual interacts with any other individual 
with equal probability. It is worth emphasizing that the 
dependence of the final state probabilities on the initial 
composition of the population is very close to the cor- 
responding quantities on finite-dimensional systems |12| . 
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This coincidence is a reflection of the conservation of the 
global magnetization by the dynamics. Indeed, for the 
classical voter model, the conservation of the magnetiza- 
tion immediately leads to the fact that the probability of 
+ consensus is equal to the initial density of + spins in 
any dimension. 

There are two extensions of the model that are worth 
mentioning. First, it is natural that extremists have a 
stronger conviction than centrists. This suggests the gen- 
eralization where a centrist adopts an extremist's opin- 
ion with probability p > 1/2 while an extremist adopts a 
centrist's opinion with probability q = l—p in an elemen- 
tary update step. For this generalization, the densities 
of each species evolve by the following rate equations in 
the mean-field limit 

x = v z x 
y = v zy 
z = —v z (x + y) 

where v = p — q quantifies the bias. For p ^ q, The 
solution to these equations are 

x(t) = m^i 

and similarly for y{t), while z(t) = 1 — x(t) — y{t). In the 
long time limit and for p > q, we find 

_ a;(0) _ y(0) 

X ^ x(0) + y(0Y y ^x(0) + y(0Y 

while z — ► 0. Thus if there is an innate bias favoring 
extremism, a frozen final state with no consensus is in- 
evitable. On the other hand, if centrists are dominant 
(p < q), there centrist consensus is the final result. 

Our results can also be adapted to determine the final 
states probabilities in the mean-field limit of the Axelrod 



model for the case of two features and two states per 
feature. The four possible states of an individual can be 
represented by 



+ + +- -+ . 

We define the dynamics to be that if an interaction pair 
has one common feature, then the other feature where 
disagreement exists is picked and the initial agent takes 
on the state of the interaction partner. 

There are four possible consensus states: all , — h, 

H — , and H — h, and two types of frozen states: mixtures of 
and H — h, and mixtures of — I- and H — . We can iden- 
tify the states and ++ in this Axelrod model with the 

— and + states, respectively, in our spin-1 Ising model, 
and the states — h and H — taken together with the 
state of the spin-1 Ising model. Then, using our previous 

results, if the initial densities of and H — h states are x 

and y respectively, the probabilities of and H — h con- 
sensus and the probability of frozen mixture of and 

H — h are given by the expressions P-(x, y), P + {x,y) and 

P_i (x,y) respectively (Eqs. I|22[) and $13$ ). By similar 

arguments, the probabilities of — h and H — consensus 
and the probability of — h and H — frozen mixture are 
given by the above expressions where x and y are now the 
initial densities of — h and H — states respectively. Un- 
fortunately, this exact mapping onto an equivalent Ising 
spin system does not work when there are more than two 
features and/or more than two states per feature. 
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APPENDIX A: SOLUTION TO EQUIVALENT SCHRODINGER EQUATION 

We present the solution to Eq. . This turns out be the Schrodinger equation for the Poschl- Teller potential hole 
|14|. for which the generic form of the equation is 



d 2 u(9) 

de 2 



2 „,2 



X(X-1) | A(l-A) 
sin 2 a9 cos 2 \9 



*(*) = 0, (Al) 



with k 2 = 1 + to 2 , a = 1, \ = A, and x(x ~ 1) = 3/4. This last equation has the roots %i = 3/2 and \2 — —1/2. The 
eigenvalues and eigenvectors of Eq. IjAlfl are (see, for example, [14] page 92): 

k 2 = a 2 ( X + \ + 2n) 2 

1 

2 ; 



u n {9) = sin x a.6 cos A a9 2F\{—n, x + A + n, x + 7S sin 2 a9), 



with n = 0, 1, 2, .., and where 2-F1 is the hypergeometric function. 

For x = ~ 1/2, u„(ff) diverges at = and = 7r/2. Thus we take the solution for x — 3/2 only and obtain the 
eigenfunctions and eigenvectors of Eq. (|llfl : 

1 + TO 2 = (3 + 2n) 2 (A2) 
u n {0) = sin^ 9 cos^ 9 2Fi(-n, 3 + n, 2; sin 2 0). 

The solutions for the angular function in Eq. Q) therefore have the form 

9„(0) = f(9) u n {9) = sin 2 29 2 F 1 (-n, 3 + n, 2; sin 2 0), 
with m given by Eq. (|A2[) . Finally, the generic form of the product solution F n (p, 9) = R n (p) O n (#) to Eq. is 

F n (p,0) = [A +P 2 ^ + A_ P - 2 ^) Q n (9). 
The coefficient A_ must be zero because F is zero at the origin. Thus 

F n {p,9) = p 2i - 2+n ^ sin 2 (2(9) 2 Fx(-n, 3 + n, 2; sin 2 (9), 
and the general solution to Eq. (JSJ) is 

00 

9) = Cn p 2{2+n) sin 2 29 2 Fi(-n, 3 + n, 2; sin 2 (9). (A3) 

Tl = 

We need the set of hypergeometric functions 2F1 to form an orthogonal set to obtain the coefficients c„ from the 
boundary condition: 

00 

F(p = 1, 9) = 1 = ^ c„ sin 2 26* 2 i<i(-n, 3 + rc, 2; sin 2 9). 

Fortunately, the functions 2-F1 are related, for certain specific integer parameters, to the Associated Legendre Poly- 
nomials (that are known to form an orthogonal set) via [l5j| : 

nm , , (-i) m r(i/ + m + i) (i-x 2 ) m / 2 _, , , i-»v 

W = 9mTV 1 2 Ji(m - 1/, m + 1/ + 1; m + 1; — - . 

2 m 1 (j/ — to + 1) to! 2 

If we take to = 1,i/ = 71 + 1 and x — cos 20 we obtain: 

— 2n' 

2 F 1 (-„, n + 3; 2; sin 2 0) = (n + 2) , gin2g (cos 20) . 
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Then the general solution Eq. (|A3|) can be written as 

oo 

F(p,8) = J2 C n Sin2d P n+l( cos2e )> (A4) 



where we absorbed the factor ( ra +2)! ^ ne coefficient c„ 



1. Solution for P+- 

To determine P^ , we now apply the boundary condition F = 1 at p = 1 and obtain 

oo 

F(p= 1,9) = 1 = E C " sin26»P,i +1 (cos26»), 

n=0 

and using the orthogonality- normalization of the P„'s we obtain the c„'s: 

n odd, 



(2n + 3) (A5) 

—7 7t n even, 

(n + 2)(n + l) 



Substituting these coefficients in Eq. (|A4(I . we obtain Eq. |1 2|l. 



2. Solution for P 4 



To determine P+, we first need the first-passage probability F(x'\x,y) for the effective random walker to hit the 
specific point (x 1 , 1 — x') on the locus x + y = 1, when starting from (x, y). This first-passage probability also satisfies 
the basic differential equation ©, but with the boundary conditions given in Eq. (|21Jl . In polar coordinates the last 
of these boundary conditions is F(x'\x, 1 — x) = F(0'\p = 1,8) = (5(cos 2 9 — cos 2 9'). Thus we now determine the 
coefficients c„ in Eq. I|A4(1 by 

00 

<5(cos 2 9 - cos 2 0') = c n sin(20) P^ +1 (cos 2(9). 

71 = 

Inverting this relation and using the orthogonality of the P^'s [Eq. l|A5(l ]. we obtain 

= (2n + 3) P n 1 +1 (co S 2gQ 
Cn (n + 2)(n + l) sin 29' 
Using this result and transforming back to xy coordinates, we obtain the first-passage probability 

v) = E ^4 4^=4, >< ^ (* + y)" ^ ( ^) • (A6) 

n(n + 1) ^/x'(l - x') \x + yj 

We then substitute this expression for F(x'\x,y) into Eq. H2()|l to determine P+(x,y). In so doing, we encounter 
the integral 

1= fx'^'-Vclx'. 

Jo y/x'il-X 1 ) 

To evaluate this integral, we first make the variable change x — 2x' — 1, use the identity Pn{x) = \/l — x 2 dP ^ x ^ , and 
integrate by parts to give 

1 {x + 1) p*(x) dx = r 1 (x + i) dp n ( x ) dx 



-1 



2 J-i 2 dx 



(x + 1) 



Pn{x) 



1 " l 



el 



-1/2 / P n {x)dx, 

-1 



Pn(l) - - y P (x)P n (x) dx = l- S n0 , (A7) 
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where we have used Pq(x) = 1, -Pn(l) = 1 Vn, as well as 

2 

P m (x) P n (x) dx = - — — 6, 



Using the results of Eq. (|A6|) and Eq. I|A7() in Eq. (|2"U|) . we obtain the solution quoted in l|22p. 

APPENDIX B: REDUCTION OF P+_ 

To reduce the series representation <|13|) for to a closed form, we start with the identity 

{0 n even, 

(-l)^n!! 
- odd < 

and make the variable change n = 2m + 1 with m > integer to obtain 

p+ _ {z) = £ (-ir(4 TO + 3)(2m -l)!! _ (m+1) 
+ w ^ 2(m+l 2m !! 1 ' y ' 

This infinite series can be summed by the following algebraic manipulation. We define w = (1 — z) 2 . The factor 
4ra + 3 can be rewritten as (2m + 1) + 2(m + 1), so that Eq. I|B1|) becomes 

p ^ = V (-l) ? "(2m + l) (2m - 1)!! (m+1) ~ (-1)" 2(m + 1) (2m - 1)!! (m+1) 
+ ~ W ^ 2(m+l)(2m)!! ^ 2(m+l)(2m)!! 

^ (2m -1)!! , - (2m-l)!! 

= ~L (2m)! , H») + W L (2m)!! (-">) > 

m=l v ' m=0 v ; 



= 1 - 



m— v 7 

1 — iu 



where we used the Taylor series expansion 

1 _ ^ (2m -1)!! 
VT=x~ ~ ^ n (2m)!! 

v m— 



This then gives closed expression for the probability of reaching the mixed state as a function of z quoted in Eq. i|14fl • 



